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Abstract. The feasibility of observing the C 18 (3-2) spectral line in cold clouds with the APEX telescope has been tested. 
As the line at 329.330 GHz lies in the wing of a strong atmospheric H2O absorption it can be observed only at high altitude 
observatories. Using the three lowest rotational levels instead of only two helps to narrow down the physical properties of dark 
clouds and globules. The centres of two C' s O maxima in the high latitude low mass star forming region CG 12 were mapped 
in C 18 (3-2) and the data were analyzed together with spectral line data from the SEST. The T M b (3-2)/T M b (2-1) ratio in 
the northern C 18 maximum, CG 12-N, is 0.8, and in the southern maximum, CG 12-S, ~2. CG 12-N is modelled as a 120" 
diameter (0.4pc) cold core with a mass of 27 M G . A small size maximum with a narrow, 0.8 kms~', C 18 (3-2) spectral line 
with a peak temperature of T M b ~ 1 1 K was detected in CG 12-S. This maximum is modelled as a 60"-80" diameter (~0.2pc) 
hot (80 K < r ex < 200 K) ~1.6 M Q clump. The source lies on the axis of a highly collimated bipolar molecular outflow near its 
driving source. This is the first detection of such a compact, warm object in a low mass star forming region. 
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1. Introduction 

The two lowest rotational transitions of 12 CO, 13 CO and 
C I8 have been used extensively for studies of dark clouds 
and globules. Also the 7 = 3-2 lines of 12 CO and 13 CO can 
be routinely observed. In contrast, there are very few pub- 
lished observations of the J - 3 - 2 transition. Besides the 
nine point C I8 (3-2) map of S106 (Little et al. 1995) and 
mapping of the centre of the pre-stellar cloud core L1689B 
( |Jessop & Ward-Thompson 200 1) only pointed observations 
are available (e.g., J0rgensen et al. 2002, 2004; Schoier et al 
(2002); Lee et al. 2003). The addition of the C ls O (3-2) tran- 
sition to the lower transitions would, however, substantially 
improve the accuracy of physical parameters derived from 
C I8 data as the relative line intensities depend heavily on the 
cloud temperature and density. 

The SEST telescope has been used to map numerous 
southern molecular clouds in C 18 (1-0) and (2-1). The new 
Atacama Pathfinder Experiment (APEX 1 ) provides an excel- 
lent opportunity to complement these data with C 18 (3-2) ob- 
servations. The APEX beam size at 329 GHz, -19", is about 
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equal to the SEST beam size at the C ls O (2-1) frequency, 24". 
This simplifies the comparison of APEX and SEST data. 

In this Letter we demonstrate the use of C ls O (3-2), (2-1) 
and (1-0) observations for radiative transfer modelling in the 
case of Cometary Globule 12, CG 12, which is a high-latitude 
star forming region at a distance of 630 pc (Williams et al. 
1977). The head of the cloud contains a highly collimated bipo- 
lar outflow (White et al. 1993) and two compact C I8 maxima 
separated by about 3' (Haikala & Olberg 2006, Paper I). The 
centres of these maxima were mapped in C 18 (3-2). Based on 
these observations a simple model is derived for both cores. 

2. Observations and data reduction 

The observations were made with APEX on Aug. 21 2005 
in good weather (PWV 0.7 mm). The C ls O (3-2) line at 
329330.546 MHz was observed with the APEX-2A SIS DSB 
receiver. The zenith optical depth was ~ 0.27 and the DSB sys- 
tem temperature ranged from 180 K to 280 K. The receiver 
signal sideband lies in the wing of an atmospheric H2O absorp- 
tion line whereas the mirror sideband centered at 335.330 GHz 
is at a more transparent frequency. The difference in the atmo- 
spheric opacities was estimated using an atmospheric model 
and was taken into account in the calibration. 

Two 5 by 5 point maps with 10" spacing were obtained 
in the position switching mode. Off positions (+5' in azimuth 
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Fig. 1. The Apex C ls O (3-2) (hanning smoothed) and SEST 
C I8 (2-1) and (1-0) (heavy histogram, light histogram and 
continuous line, respectively) spectra in CG 12-N. The map 
grid step is 10". The tick mark is at velocity -5.9 kms -1 and 
the map centre position is 13 h 57 m 39 s !5, -39°56'3" (J2000). 

for CG 12-N and +5' in azimuth and +5' in elevation for 
CG 12-S) were outside the globule at the time of the obser- 
vations. CG 12-N was observed at elevations 60° to 35° and 
CG 12-S at 32° to 22°. Integration time of 20s was used and 
calibration was done every 10 to 30 minutes. Each position was 
observed two or three times. The map centre positions were in- 
tegrated for a total of 250s. The 1 GHz bandwidth of the MPIfR 
Fourier transform spectrometer was divided into 16384 chan- 
nels of 61 kHz ( -55ms 1 at 329 GHz). 

Most of the observed spectra contain ripple due to varia- 
tions in the atmospheric emission, reflections in the telescope 
optics and instability of the receiver. In the data reduction the 
possible low frequency ripple was first fit with a sinusoidal 
baseline where-after possible higher frequency ripple was re- 
moved by masking the corresponding frequency in Fourier 
transform space if possible. Finally a first order baseline was fit 
to the spectra around the source velocity and subtracted. The re- 
sulting rms of the spectra is by ~15% higher than the value ex- 
pected according to the radiometer formula indicating that the 
possible non white noise due to, e.g., remaining baseline ripple, 
does not dominate the spectral noise. The spectra were trans- 
formed to the main beam brightness temperature scale using 
a main beam efficiency of 0.7. All the spectral line intensities 
and line integrals reported in this paper are on the Tmb scale. 

3. Results 

The hanning smoothed C ls O (3-2) spectra in CG 12-N and 
CG 12-S are shown in Figs. ^ an d 12 together with the SEST 
C 18 (2-1) and (1-0) observations (Paper 1). The area mapped 
in C 18 (3-2) in each maximum covers only the SEST HPBW 
at 109 GHz, ~47", and the emission continues beyond the map 
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Fig. 2. As Fig.[T]for CG 12-S. The tick line is at -6.2 kms" 1 
and the map centre position is 13 h 57 ra 41 s .4, -39°59'3" (J2000). 
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Fig. 4. C 18 , CS (3-2) and DCO + spectra in the C 18 (3-2) 
map centre positions. For CG 12-S also a CS(5-4) spectrum is 
shown 

boundaries. Maps of the C ls O (3-2) line integral are shown in 
Fig EJsuperposed on SEST molecular line and mm continuum 
maps and a NIR Ks (NTT SOFI) image. 

CG12-N: In the southern part of the CG 12-N map the 
line profiles and intensities of all the observed C I8 transitions 
are similar. Elsewhere both the C ls O (3-2) line profile and 
line peak velocity vary more than in the (2-1) transition. The 
maximum is shifted to SE from the C I8 (2-1) (and also (1- 
0)) maximum and is located near the faint IRAS point source 
12546-3941 (Fig.|3}. 

The three C ls O transitions together with the CS (3-2) and 
DCO + (2-1) lines in the CG 12-N map centre position are 
shown in Fig. [4] The C I8 lines are asymmetric but the pro- 
files are similar. An unpublished C I7 (1-0) spectrum obtained 
in this position shows that the C ls O (1-0) line is optically thin. 
The weak (Tmb ~ 0.5K) CS line peaks at the same velocity as 
the C 18 lines. The DCO + line is redshifted with respect to the 
C ls O and CS lines. 

CG12-S: The C I8 (3-2) emission peaks 10" West of the 
CG 12-S map centre position. In contrast with CG 12-N, the 
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Fig. 3. C 18 (3-2) line integral and SEST mm line and continuum (arbitrary contour levels) and NTT NIR Ks observations. 
The upper row shows CG 12-N and lower row CG 12-S. C ls O (3-2) is shown in gray scale except in the NIR Ks panels where 
contours are used. Lowest C I8 (3-2) contour level and the increment are 2.0 Kkms 1 and 0.6 Kkms 1 for CG 12-N and 3.8 
Kkm s _1 and 0.9 Kkm s _1 for CG 12-S. The arrow in the CG 12-S Ks panel shows the orientation of the outflow axis and points 
at the assumed position of the driving source. The APEX and SEST beam sizes are shown in the upper left corner of each panel. 
The offset from the map centres in Figs.^and|2]in arc seconds is shown on the axes. In the rightmost panel the SEST C 18 (2-1) 
emission contours superposed on the SOFI Ks image is shown. The lowest contour is 1.6 Kkms 1 and the increment is 0.6 
Kkms~'. The origin is that of Fig.[2]and the box shows the location and size of the CG 12-N C 18 (3-2) map. The red and 
blue outflow lobes from White (1993) are shown with heavy black and white contours, respectively. The positional uncertainty 
ellipses of the IRAS point sources 13547-3944 (in white) and 13546-3941 (in black) are also shown. 



C ls O (3-2) line is significantly stronger than C ls O (2-1) in the 
map centre, West and South- West. The maximum lies at the 
NW edge of the dense core outlined by the high density tracers 
(DCO + (2-1), CS (3-2)) and the 1.2 mm continuum and on the 
axis of the outflow. It is likely that the driving source of the 
molecular outflow (White 1993) lies at the Southern tip of the 
cone-like NIR nebulosity seen in the Ks image at +5", -20". 
The nature of the driving source is not known. 

In the C ls O (3-2) map centre position two C I8 velocity 
components (-6.4 km s~' and -6.2 km s -1 , Paper 1) are blended 
and the lines are skewed (Fig.|4jl. Only the -6.4 kms~' com- 
ponent emits strongly in DCO + , CS (2-1), CS (3-2) and 
H 13 CO + (1-0) (Paper 1). However, in C 18 Othe -6.2 kms 1 
component is the strongest. A two component gaussian fit with 
centre velocities fixed to -6.4 kms~' and -6.2 kms~' gives 
7.6 K, 3.7 K and 1.7 K, respectively, for the C ls O (3-2), (2- 
1) and (1-0) peak temperatures for the component centred at 
-6.2 km s _1 . In C 18 (3-2) this component extends beyond the 
mapped area in the South- West (Fig|2j. The source is elongated 
in the NE-SW direction but is hardly resolved perpendicular to 
this direction. 

4. Discussion 

The T MB (3-2)/T MB (2-1) ratio (R 32 ) is 0.8 and 2 (the -6.2 
kms 1 component) in CG 12-N and CG 12-S , respectively. 
Assuming LTE, these ratios would correspond to T ex of <10 K 
in CG 12-N and >30 K in CG 12-S .Many of the C 18 Q line ra- 
tios in pre- and protostellar cores ( J0rgen sen et al. (2005)| and 



references therein) are similar to those observed in CG 12-N. 
These line ratios were modelled with a "drop" model which 
assumes a strong molecular depletion in a cold shell (<30 K) 
between the heating source and the outermost part of the en- 
velope. The intensity of the C 18 lines in IRAS 16293-2422 
is as high as observed in CG 12-S and R32 is 1.3 (Schoier 
et al, 2003). The lines are, however, optically thick and, 
e.g., the CS(5-4) line is stronger than the C 18 (3-2) line. In 
CG 12-S the CS(5-4) line is not detected (Fig.gJ. 

It is evident from Figs.^and[2]that considerable fine struc- 
ture in C 18 (3-2) is present both in intensity and velocity in 
the two mapped areas. However, the signal-to-noise ratio of the 
spectra is not high enough for a more refined analysis of the 
source structure than presented in Sect. [5] We therefore restrict 
our further analysis to constructing models which reproduce 
the observed line ratios and intensities in the C ls O (3-2) map 
central positions (Fig.|4ji where the signal-to-noise ratio is best. 
In CG 12-S only the -6.2 kms 1 component is considered. 

The reproducibility of the C 18 (1-0) and (2-1) spectral 
line intensities at SEST has been very good. The lines were 
also observed simultaneously and therefore the observed rela- 
tive line intensities should be quite accurate. No such record is 
available for the APEX telescope and the receiver. We therefore 
assume a calibration accuracy of 20%. 

Radiative transfer modelling was done with the Monte 
Carlo method (Juvela 1997) assuming a spherically symmet- 
ric geometry. The C 18 maps were used to estimate the cloud 
radii, 60" for CG 12-N and 35" for CG 12-S. The density dis- 
tribution adopted as a starting point was n oc r~ L5 truncated at 
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cloud radius. Power law temperature distributions were tested 
with temperature either increasing or decreasing towards the 
cloud centre. The model was fitted to the intensities and line 
widths of the observed C 18 spectra (CG 12-N) and to line 
temperatures fitted for the -6.2 km s -1 component (CG 12-S ). 
The spectra around the central position were used to constrain 
the radial density profile. The free parameters of the fit were 
the size of the inner region with constant density and temper- 
ature, and the parameters of linear scalings that were applied 
separately to the density and temperature values. The model 
spectra were convolved with gaussian beams that correspond to 
the resolution of the observations. Constant C 18 abundance of 
2.5 10 -7 was used. The quality of the fit depends only weakly 
on the abundance which does, on the other hand, directly affect 
the mass estimates. 

The resulting radial density and temperature profiles are 
shown in Fig. 5. For CG 12-N the ^ 2 -value is close to one. 
However, in the case of CG 12-S the ;^ 2 -value is over five, 
mainly because the predicted 7=1-0 intensity is only 
half of the observed value. Best fits in CG 12-S were obtained 
with models where the kinetic temperature increased inwards. 
One could obtain reasonable fits even with isothermal mod- 
els. CG 12-N is modelled with a 120" diameter cloud where 
the temperature decreases from 12 K down to 6 K at the edge. 
The CG 12-S data can only be fit with a small size, hot clump. 
The modelled total gas column density in the APEX beam and 
the cloud mass are 4.010 22 cirT 2 and 27 M Q for CG 12-N and 
7.810 21 cnr 2 and 1.6 M Q for CG 12-S, respectively. 

The models presented above should be considered only 
as indicative. Spherical symmetry was assumed and the 
C 18 (3-2) maps cover only the very centres of the maxima. 
Furthermore, the cloud sizes are close to the telescope beam 
size so that the beam convolution has a large effect on observed 
and calculated line ratios. Especially the CG 12-S fit is sensi- 
tive to the beam convolution. Different combinations of source 
size, temperature and density produce acceptable fits. The cen- 
tral density can not be increased much as the CS lines would 
become stronger than observed. Common to all of these solu- 
tions is the small size, a radius of 30" to 40", and high cen- 
tral temperature, 80 K to 200 K. The problem is, however, not 
unique for this source but inherent to all observing and mod- 
eling small size (compared to beam size) sources. The source 
beamfilling factor becomes an important modelling parameter. 

If the outer parts of the CG 12-N cloud envelope are heated 
up, e.g, by interstellar radiation field (e.g., Snell 1981), this has 
very small effect to the fit at the cloud centre position. The gas 
column density in the cloud envelope is small compared to the 
central density from where the bulk of the emission is coming. 
The LTE mass ( r ex = 10 K) calculated from the C 18 (2-1) 
data within 60" from the centre of CG 12-N is 25 M which is 
in agreement with the modeled mass of 27 M . 

The LTE mass of the -6.2 km s _1 velocity component in 
CG 12-S calculated from C 18 (2-1) data within 35" from the 
C 18 (3-2) map centre is 4.1 M Q for r ex =15 K. However, for 
r ex =160 K the mass would be three times higher. The LTE 
method assumes that all energy levels are thermalized and this 
is clearly not the case in the CG 12-S model. 
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Fig. 5. The radial density (lower curve) and temperature pro- 
files of CG 12-N and CG 12-S. The horizontal axis shows the 
distance from the clump centre. One arc-min corresponds to 
0.18 pc or 3.8 10 4 AU at the assumed distance of 630 pc. 

5. Conclusions 

The addition of the C 18 (3-2) transition to the C ls O data on 
the two lower transitions restricts strongly the parameter space 
of possible cloud models. 

The C 18 (3-2) observations of CG 12-N confirm the cold 
temperature of this core. The C ls O (3-2) line profile varies 
more over the mapped region than the C ls O (2-1) profile. The 
C 18 (3-2) line is therefore probably better suited for study- 
ing the small-scale stucture of molecular clouds than the lower 
transitions of this molecule. 

The detection of high, -10 K, C 18 (3-2) line tempera- 
tures in CG 12-S was unexpected and forced us to reconsider 
the temperature and density structure of the material traced 
by C 18 (2-1) and (1-0) presented in paper 1. The small, hot 
clump lies projected on the edge of a dense core traced by 
DCO + . It also lies on the axis of a highly collimated bipolar 
molecular outflow near its driving source. This suggests that the 
two are related. As the C ls O lines from the hot spot are narrow 
and appear at the same radial velocity as the parent cloud, they 
probably do not originate in shocked gas. Neither is direct ra- 
diative heating from the strong IRAS point source 13547-3944 
likely. The heating mechanism and the relation of the hot spot 
to the collimated outflow remain therefore unspecified. Further 
studies of this phenomenom seem warranted as they may lead 
to a better understanding of the interaction between newly born 
stars and the surrounding molecular material. 

Better spatial resolution can be achieved by observing 
still higher C 18 transitions which should be detectable in 
CG 12-S. The sizes of the CG 12 C ls O maxima are well suited 
for the future ALMA interferometer 
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